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^ ! We study the possibihty of the scalar charmonium and glueball mixing in 

I e~^e~ annihilation at ^/s = 10.6 GeV. The effects can be used to explain the 

^ I unexpected large cross section (12 ± 4 fb) and the anomalous angular distri- 

■~-«^ ■ bution (a = — I.IIqJ) of the exclusive e~^e^ J/'ipXco process observed by 

• Belle experiments at KEKB. We calculate the helicity amplitudes for the process 

ji^! e+e" — > J/tljH{0^~^) in NRQCD, where i7(0++) is the mixed state. We present 

D I a detailed analysis on the total cross section and various angular asymmetries 

' which could be useful to reveal the existence of the scalar glueball state. 

> ■ 
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1 Introduction 

The charmonium production has long been considered as a good process for investigating 
both perturbative and nonperturbative properties of quantum chromodynamics (QCD), be- 
cause of the relatively large difference between the scale at which the charm-quark pair is 
produced at the parton level and the scale at which it evolves into a quarkonium. In partic- 
ular, comparing to hadron colliders, e^e~ colliders, provide a cleaner environment to study 
the charmonium productions and decays. However, some puzzles arise from the recent mea- 
surements on the prompt productions at BaBar and Belle [1, 2, 3]. For the inclusive J/ if) 
productions, the cross section is much larger than the predictions of nonrelativistic quantum 
chromodynamics (NRQCD) [4]; there is also an over- abundance of the four-charm-quark 
processes including the exclusive J/ip and charmonium productions; there is no apparent 
signal in the hard J/ip spectrum which has been predicted by the J/ipgg production mode 
as well as the color-octet mechanism in NRQCD. To provide plausible solutions and expla- 
nations for these conflicts, theorists have studied the possibilities of the contribution from 
two- virtual-photon mediate processes [5], large higher-order QCD corrections [6, 7], collinear 
suppression at the end-point region of the J/ip momentum [8, 7], contribution from the 
J/-?/'-glueball associated production [9] and contribution from a very light scalar boson [10]. 

Although all of the above corrections significantly enhance the cross sections for the 
prompt J/ip productions, they can be well distinguished by their momentum and angular 
distributions. The two- virtual-photon processes have a very hard spectrum for the J/ijj 
momentum and an enhancement in the large |cos^^| region (6 is the J/ip scattering angle 
in the ccnter-of-mass frame) due to the t-channcl electron-exchange contribution [5]. The 
collinear suppression softens the hard spectrum of the color-octet production [8] as well 
as that of the color-singlet J/i/^gg process [7], while the high-order QCD corrections are 
generally expected not to change the angular distribution strongly and their effect could be 
represented by a large renormalization K factor [7]. 

J/-?/'-glueball associated production in e+e^ annihilation has been shown to be of a spe- 
cial interest recently [9]. As a basic prediction of QCD, glueball states have been studied 
in the QCD sum rules [11, 12], bag model [13, 14], nonrelativistic potential model [15], and 
lattice QCD [16]. In the search of the glueballs, radiative decays of heavy quarkonium, such 
as J/ip '-fX and T 7X, are regarded as favorable processes for glueball production 
and high statistics data have been studied by BES and CLEO Collaborations [17, 18]. The 
quarkonium-glueball associated productions at B factories should be another good place to 
discover glueballs because their amplitudes are essentially the same as those of the quarko- 
nium decays. If glueball masses lie in the charmonium mass region, the glueball may be 
hidden in the charmonium resonance peak in the exclusive production, and this mechanism 
may contribute to the excess of exclusive charmonium-pair production events. In Ref. [9], 
J/'i/'-glueball productions e~^e~ — > J/ipgg J/ipQ have been studied, and at the leading- 
twist level and the leading-order of the NRQCD velocity-power-counting rule, the production 
of the scalar glueball state ^o(0"'"''") is found to be dominant. As we shall see, the angular 
distribution of e'^e' — > J/ipQo is dramatically different from those of the other exclusive 
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J/'0-charmonium productions, such as e'^e~ — > J/iprjc and e'^e~ — > J/ipXco- The angular 
distributions of all those processes with single virtual photon exchange can be parametrized 
as 

+ 1 

a cost* 

by using an asymmetry parameter a (—1 < « < 1). For e~^e~ J/ipfJc and e^e~ J/ipXco, 
a = 1 and 0.25 respectively [19], while for the J/ip-Qo production mode, a is estimated to be 
—0.85 in NRQCD [9]. We note here that preliminary result on the J/ijj angular distribution 
in the J/ip-Xco region is o; = — Lllgg, which is more consistent with o; ~ — 1 rather than 
a ~ 0.25 [3]. Conservatively, the present Belle data [3] suggest that a should be negative. 
This unexpected angular distribution leads us to consider further the possibility of the scalar 
charmonium Xco and glueball mixing. In this Letter, we will focus on the mixing effects 
on various angular distributions which could be useful to reveal the existence of the scalar 
elucball state. 



2 Formulae 

Let us start with the angular momentum and parity considerations. The J/ilj-scalai{H) 
associated production via a virtual photon in the J-^*^ notation is 1 1 + 0++. If we 

denote by L the orbital angular momentum between J/tjj and H, the angular momentum 
conservation tells L = 0, 1, 2. Conservation of parity gives —1 = lx(— l)x(— 1)^, and hence 
L should be even. We are left with L = or 2, that is, only S-wme and D-wme J/t/^-H are 
allowed. 

Now we present the hclicity amplitudes for the J/ip-Xco and J/ip-Qo productions sepa- 
rately in NRQCD. One of the crossed Feynman diagrams for e+e^ — > J/ipXco is shown in 
Fig. 1(a), and that for e+e" — > J/ipQo is shown in Fig. 1(b). We adopt the standard co- 
variant projection formalism for the color singlet productions [20] . If a charmonium carries 
the momentum P, then the charm quark and anti-quark within the charmonium are on-shell 
and carry momenta P/2 — p and P/2 + p, where p is the relative momentum of the quark 
pair and satisfies p ■ P = 0. The amplitudes for the charm-pair production are furthermore 
projected into a certain spin component by the projection operator 

= E ^(^/2 - P; S2)u{P/2+p; si)(l/2, s^■, 1/2, s^lS, S,). (2) 

sl,s2 

After expanding the amplitudes in p and convoluting them with the wave functions at the 
origin, one obtains the helicity amplitudes for S'-wave and P-wave charmonia, such as J/ip 
and Xco- For the glueball production as shown in Fig. 1(b), we follow the method of 
Ref. [9] and extract the leading-twist contribution by expanding the gluon momenta in the 
light-cone direction and requiring the gluon pair to be collinear. We define the light-cone 
direction as = (1,0,0, 1), while the J/ip momentum direction is along n = (1,0,0, — 1). 
The formation of the scalar glueball Qo from the gluon pair is described by a collinear wave 
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(«) {b) 

Figure 1: (a) One of the four Feynman diagrams for e~^e~ — >■ J/i/jXco] (b) One of the six 
Feynman diagrams for e~^e~ J/ipQo- 

function Jq = /q^ dx(j)Q{x), where 0o is the distribution wave function for Qq and x is the hght- 
cone momentum fraction. We denote the Qq momentum by /c = {k^ = k ■ n, k~ = k ■ n, 0±), 
and further represent the transverse momentum of a gluon inside Qq as k^. By integrating 
over and transforming from the momentum to the coordinate space, the distribution wave 
function 0o(x) is obtained as 

Here = [—ga/3 + + narip) /2] / \/2 is the scalar projector. G~ is the gluon field along 
the light-cone direction and a is the color index. 

Under the above approximations, the helicity amplitudes for J/ip-Xco production are 
given by, 

16 r 1 
M^(a, A) = gyq,^^—- (Os^ - hQsmDj,- e* + (-14s + 2Aml)3,- P,e*y $(0)$'(0), (4) 

and those for J/ip-Qo production are 

M,(a, A) = gyq}^^ [-^ - ^-^ifAl^] $(0)/o, (5) 

s — Am^ 2s J 

where j is the electron-positron current and e is the polarization vector of J/ip. o and A 
denote the helicities of the electron-positron current (along the electron beam direction) and 
J/?/', respectively. $(0) and $'(0) are the wave functions at the origin for J j^l) and Xco- In 
the center-of-mass frame, the amplitudes are reduced to 

M^,g((T = ±, A = 0) = gl(?Qca^^g sin d, 



M,,,(a = ±,A = ±) = gyQJo J^^jL^'\ (6) 



where the coefficients a^^g and h^^g are for the longitudinally and transversely polarized J/ 
amplitudes respectively. We find 

32 1 



«x = C^F^^(4 + 10r2-3r^)$(0)$'(0), 



3 



S2 1 

^ ^^F^^(-18 + 7r2)$(0)<|.'(0), (7) 



for J/i^-Xdd production, and 



I S2 1 

6, = -^^-^$(0)/o, (8) 



for J/ip-Qo production. Here r = irric/ y/s and = (iV^ — l)/(2A^c) = 4/3 is the color 
factor. These results are consistent with those in Refs. [9, 19, 21]. 

It is natural to expect the mixing between the pure states of Xdo ^-nd Qq if they carry the 
same quantum numbers J^*^ = 0''"''" and have small mass difference. The pure charmonium 
and glueball states, Xco and Qq, and the mass eigenstates Hi, H2 are then linked by 

r")=^(^0' f^=('°'^. (9) 
Go J \H2J V-sm^ cos^y' ^' 

where ^ is the mixing angle. The mass-squared matrix for Xco and Qq has the following form 
and is diagonalized as 

,({M^f A \^_({M,f \ 
^ V A {M,f)^-\ {M,f)- ^^"^ 

Wc assume that the observed Xco state is Hi, and hence Mi = M(xco) observed = 3.415 GeV [22]. 
Because Hi properties observed in charmonium radiative decays are consistent with the pure 
Xco assumption, we expect the mixing angle to be small sin^^ < 0.1. The mixing angle ^ is 
real if the transition matrix element A is real. If > Mg, Mi > M2. while if M^ < Mg, 
Ml < M2. The mixing angle can be non-negligible (sin^,^ ~ 0.1) if |A| ~ — M||. 

The helicity amplitudes for J/ip-Hi and J/1P-H2 productions keep the forms in Eq. (6) 
and the coefficients a^^g and b^^g have to be replaced by aHi,H2 and bHi,H2 respectively. With 
Eq. (9), we have 

ajy^ = cos,^ — sin^, a/^j = sin^ + a^, cos^, .^^^ 

bni = &x "^^^ ^ ~ bg sin ^, bffi = b^ sin ^ + bg cos ^. 

In the sin^ — limit. Hi is a pure Xco, while H2 is a pure glueball. 

3 Results 

For J/ijj productions in e^e~ annihilation, the mixing effects can be measured via the angular 
distributions of J/ip and its leptonic decays. According to the helicity amplitudes in Eq. 
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(6), the triple angular distribution for the process e+e — > J/ipH — > I'^l H {H denotes Hi 
or H2) is given by 

da SgyQlVT^B 

- X 



, . o - |6/^P(l + cos2^)(l + cos2r) +2|aj^Psin2^sin2r 

d cos ed cos e*d(f)* 20 487r2s Li ^1 ^ ^ i 

+ 2Re(a^6j^)sin2^sin2rcos0* + |6//|^sin2^sin2rcos2(/)*], (12) 

where 6* and 0* are the polar and azimuthal angles of /~ in the J/ip rest frame. The polar 
angle 6* is measured from the J/iJj momentum direction in the e~^e~ collision rest frame, 
and (f)* is measured from the scattering plane. In Eq. (12), B is the branching fraction 
B{J/7j; in-) = 5.9% for l^e,n. 

With respect to the mass difference of the mass eigenstates Hi and H2-, we meet with 
two different cases: (1) If the mass difference is large enough, only Hi has been identified 
as the Xco resonance. (2) It is also possible that the mass difference is quite small and both 
Hi and H2 contribute to the observed resonance peak. Therefore we will consider these two 
cases separately. 

In the first case, we rename Hi as H. By integrating over 9* and 0* and comparing with 
Eq. (1), we obtain 

b%-a% aT - 2a L , . 

where aT and (Jl are the cross sections for the transversely and longitudinally polarized J/ip 
productions, and a — ar+aL is the total cross section. One can obtain the cos 9* distribution 

by integrating over 6 and </>* in Eq. (12), and obtain da/dcos9* ~ 1 +a*cos'^6*. It is easy 
to find a* = a according to Eq. (12). The interference term 2Re(a^6j/) can be measured 
through the combination of the 9, 9* and (p* distributions. We introduce a new quantity 
as the normalized interference term 

^ ^n^ ■ ^^^^ 

Both a and arc bounded in the interval [—1, 1]. 

Now we present our numerical results based on the input parameter values e^/47r = 
1/129.6, qI/At: = 0.26, = 1.5 GeV, and = 10.58 GeV. The value of the J/-^ wave 
function at the origin can be obtained from V{J/^p e+e'), |$(0)p = 0.0336 GeV^ in the 
leading order of NRQCD [19, 6, 7]. 

The Xco wave function at the origin $'(0) is estimated from the observed width T{xco 
77) [19]. Since we assume that the observed state is a mixture H, the partial width is 
actually proportional to [ cos^$'(0)|^. Here we assume that only the charmonium content 
of H contributes to the decay Xco 77 at the leading order in 0:3. Therefore we obtain 
I cos ,^$'(0)1^ = 0.0117 GeV^. It should be noted that the partial widths and the observed 
radiative transition rates can be affected at several percent level for the mixing angle of the 
order of sin^ ^ ~ 0.1. 
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The upper bound of the glueball wave function |7o| is obtained from the radiative T decay. 
The CUSB Collaboration reported a 90%-C.L upper limit for the branching ratio Br{T 
^H) which is about 0.01%-0.15% for the scalar boson (H) mass between 2 and 8.5 GeV [23]. 
The upper bound is obtained in terms of the CUSB mass resolution, 20 MeV. As discussed 
in Ref . [9] , if the H decay width Th is larger than the resolution, the upper bound should 
be loosen by a factor of r7j/(20MeV). In Ref. [9], the authors considered the pure glueball 
production process, e~^e~ — >■ J/ijjQo, and they assumed the glueball decay width should be 
less than about 100 MeV, which is the full width at half maximum of the 'x^o' peak in the 
Belle fit to the J/tp recoil mass distribution. They found the upper bound for the glueball 
wave function |/op < 5.8 x 10~^ GeV^. Accordingly, the rate of the cross sections of the pure 
glueball and charmonium productions crj/^g^/aj/^^^Q which is proportional to |/oP/|^'(0)P, 
is less than 0.72. Similar upper bound applies in our case, except that the radiative T decay 
rate constrains | sin,^/op instead of |/oP, since T can decay to the glueball content of H only, 
at the leading order in a^. In the e~^e~ J/ipH process, because the contributions from 
the charmonium and glueball contents are proportional to cos^ ^ and sin^ ^, respectively, 
the upper bound of Ref. [9] leads to the bound (sin^^(Tj/^g(,)/(cos^^(7j/^j^^o) < 0.72. We 
therefore define the following quantity 



7^ = Sign 



sin^- 



$'(0)J COS^e (7j/v,x.o 



to present our numerical results. Here Sign 



lo 



determines the sign of the interference 



terms. The ^-dependence of our results enters only through the rate TZ. For simplicity, we 
assume that the wave functions $'(0) and Iq are real. 

The magnitude of the parameter TZ can be as large as unity if Hi is the only gluon-rich 
state which contribute to the radiative T decay. It is, however, unlikely that the H2 mass 
lies outside of the reach of the CUSB search [23]. If the H2 mass is in the region 2-8.5 GeV, 
the radiative T decay rate should constrain | cos^/op, and therefore gives an upper bound 
on |7^| in Eq. (15) of the order of tan^ CirH^/WO MeV) x 0.72. If sin^ C < 0-1 and if the H2 
width r^/j is below 500 MeV, the allowed range of TZ is reduced to \TZ\ < 0.4. 

In Fig. 2, we show the cross section for e^e" — * J/ipH in the interval —1 < TZ < 1. We 
find the cross section from the contribution of the Xco content, cos^,^ '^J/iPxca^ 2.7 fl:), much 
smaller than 12 ±4 fb which we estimate from the preliminary Belle data [24]. We note that 
the contribution from the interference term can significantly increase the cross section in the 
region — 1 < 7?. < 0, and decrease the cross section for < 7^ < 1. For instance, at 7^ = —0.1, 
the total cross section is about 4.3 fb, while sin^^ ^J/'ipQa = 1"^! cos^ ^J/ipxco — 0-27 fb. The 
interference term makes up the difference 4.3 — 2.7 — 0.27 = 1.3 fb of the total cross section, 
which is about five times greater than the direct contribution from the glueball content. 
At TZ — —0.4, the total cross section is about 6.5 fb among which 2.7 fb comes from the 
interference contribution. While at 7^ = — 1, the total cross section can be as large as 9.6 
fb (4.2 fb from the interference contribution), which is comparable with the central value of 
the experimental data. 

In Fig. 3, we plot the angular asymmetry parameters a (solid line) and ctoff (dashed line) 
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Figure 2: The cross section for e~^e~ — * J/ipH in the interval —1 < TZ < 1 for the case that 
only one of the mass eigenstates lies in the 'Xco' resonance. 
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Figure 3: The asymmetries a and a^s in the interval —1<TZ<1 for the case that only one 
of the mass eigenstates lies in the 'Xco' resonance. 

which are defined in Eqs. (13) and (14) as functions of TZ. a for the charmonium content 
is 0.25 at 7^ = 0, while for the glueball content is -0.85. The Belle fit finds a = -l.llae, 
consistent with —1 [3]. Although the error is still quite large, we may state that a < is 
favored by the Belle data [3]. In Fig. 3, « < corresponds to the region TZ < —0.045 and 
TZ > 0.85. At 7^ = —0.4 and TZ = —1, a can be as small as —0.30 and —0.43, respectively. 
The region TZ > 0.85 is not favored because of the suppressed cross section in Fig. 2. 
For the case of the asymmetry parameter aos, (Xos — —0.97 for the charmonium content 
{TZ = 0), and aos — —0.52 for the glueball content {TZ — > ±oo). According to Eqs. (13) 
and (14), a^s = ±1 should correspond to a = which reflects the produced J/ip mesons 
are unpolarized (ctt = '2o'l)- (^oS = corresponds to a = 1 or —1, that is, all J/ip are 
transversely or longitudinally polarized. These can be observed in Fig. 3 at 7^ = —0.045 
{a — while ctoff = —1), TZ — 0.27 {a — 1 while ctofr = 0), and TZ — 0.85 {a — while 
CKoff = 1)- In the region —1<TZ< —0.045 favored by the Belle data, ctofr varies slightly from 
-1 to -0.90. 
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Figure 4: The asymmetries a. and aofr in the interval 7^ = <^.j/i!go/ <yj/i}XcQ < ^ 
that both of the mass eigenstates Hi and H2 lie in the 'Xco' resonance. 



Now we turn to the second case in which both of the two mass eigenstates Hi and H2 are 
hidden in the 'xco' resonance peak of the J/tl) recoil mass distribution in the Belle data [3]. 
In this case, one can obtain the total cross section and asymmetries by summing over the 
triple angular distributions in Eq. (12) for Hi and H2- By using the relation between the 
coefficients aHi,H2i^Hi,H2 ^^^d a'x,gi^x,9 ^1- (1^)' ^^^^ ^^^^ 



(16) 



From the above relations, we can find that the sum of the contributions from the two mixed 
states Hi and H2 is essentially identical to the sum of the contributions from the pure Xco 
and Qo states. 

The sum of the cross sections for J/ip-Xco a-nd J/ip-Qo productions is simply to be + 
(^J/fSo = (1 + ^)<^J/V'Xco> where 



7^ = 



(17) 



For the leading order estimate of crj/^y^^^ — 2.7 fb, the sum of the cross sections can be as 
large as 5.4 fb for 1Z = 1, which is to be compared with our estimate of 12 ± 4 fb [24]. 
The asymmetry a for this case is given by 



a 



bl 



and CKofr is expressed as 



bl + al + bl + al' 



2Re{a^b^ + alb. 



(18) 



bl + al + bl + al^ 



(19) 
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where a^^g and b-^^g follow from Eqs. (7) and (8). They are plotted against 1Z in Fig. 4. 
One can find a — 0.25 and = —0.97 by setting = for the J/ip- Xco production, and 
a = —0.85 and aofr = —0.52 for the J/ip-Qo production in the TZ ±oo limit, which are the 
same as those for the previous case, a falls off and aoff goes up with increasing TZ. TZ > 0.20 
gives a < and ctoff > -0.87. At = 1, a = -0.41 and ctoff = -0.70. 

In this paper, we studied the disagreement between the experiment and the NRQCD 
prediction for the J/ip-Xco associated production, in the production cross section and the J/ip 
angular distribution. We introduced a charmonium-glueball mixing mechanism to explain 
this discrepancy and propose to study various angular distributions to explore the possible 
mixing effects. Our results may facilitate the present and future experimental measurements 
to resolve possible glueball contents in the charmonium resonances. 
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